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Abstract. The Jeans equations relate the second-order velocity moments to the den- 
sity and potential of a stellar system. For general three-dimensional stellar systems, 
there are three equations, but these are not very helpful, as they contain six inde- 
pendent moments. By assuming that the potential is triaxial and of separable Stackel 
form, the mixed moments vanish in confocal ellipsoidal coordinates. The three Jeans 
equations and three remaining non- vanishing moments form a closed system of three 
highly-symmetric coupled first-order partial differential equations in three variables. 
They were first derived by Lynden-Bell in 1960, but have resisted solution by standard 
methods. Here we present the general solution by superposition of singular solutions. 



1 Introduction 

Much has been learned about the mass distribution and internal dynamics of 
galaxies by modeling their observed kinematics with solutions of the Jeans equa- 
tions (e.g. [4]). The Jeans equations connect the second-order velocity moments 
(or the velocity dispersions, if the mean streaming motion is known) directly to 
the density and the gravitational potential of the galaxy, without the need to 
know the phase-space distribution function /. In nearly all cases there are fewer 
Jeans equations than velocity moments, so that additional assumptions have 
to be made about the degree of anisotropy. Furthermore, the resulting second 
moments may not correspond to a physical distribution function / > 0. These 
significant drawbacks have not prevented wide application of the Jeans approach 
to the kinematics of spherical and axisymmetric galaxies. Many (components of) 
galaxies have triaxial shapes ([2], [3]), including early- type bulges, bars, and gi- 
ant elliptical galaxies. In this geometry, there are three Jeans equations, but little 
use has been made of them, as they contain six independent second moments, 
three of which have to be chosen ad- hoc (see e.g. [11]). 

An exception is provided by the special set of triaxial mass models that have 
a gravitational potential of Stackel form. In these systems, the Hamilton-Jacobi 
equation separates in confocal ellipsoidal coordinates ([19]), so that all orbits 
have three exact integrals of motion, which are quadratic in the velocities. The 
three mixed second-order velocity moments vanish, so that the three Jeans equa- 
tions for the three remaining second moments form a closed system. Lynden-Bell 
([13]) was the first to derive the explicit form of these Jeans equations. He showed 
that they constitute a highly symmetric set of three first-order partial differen- 
tial equations for three unknowns, each of which is a function of the ellipsoidal 
coordinates, but he did not derive solutions. 
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When it was realized that the orbital structure in the triaxial Stackel models 

is very similar to that in numerical models for triaxial galaxies with cores ([6], 
[16]), interest in the second moments increased, and the Jeans equations were 
solved for a number of special cases. These include the axisymmetric limits 
and elliptic discs ([8], [10]), triaxial galaxies with only thin tube orbits ([12]), 
and the scale-free limit ([11]). In all these cases the equations simplify to a 
two-dimensional problem, which can be solved with standard techniques after 
transforming two first-order equations into a single second-order equation in 
one dependent variable. However, these techniques do not carry over to a single 
third-order equation in one dependent variable, which is the best that one could 
expect to have in the general case. As a result, the latter has remained unsolved. 

We have solved the two-dimensional case with an alternative solution method, 
which does not use the standard approach, but instead uses superposition of 
singular solutions. This approach can be extended to three dimensions, and 
provides the general solution for the triaxial case in closed form. We present 
the detailed solution method elsewhere ([23]), and here we summarise the main 
results. In ongoing work we will apply our solutions, and will use them together 
with the mean streaming motions ([20]) to study the properties of the observed 
velocity and dispersion fields of triaxial galaxies. 



2 The Jeans equations for separable models 

We define confocal ellipsoidal coordinates (A, /x, v) as the three roots for r of 



T + a T + f3 T + 7 

with (x, y, z) the usual Cartesian coordinates, and with constants a, j3 and 7 such 
that —7 < u < —(3 < n < —a < A. Surfaces of constant A are ellipsoids, and 
surfaces of constant ^ and v are hyperboloids of one and two sheets, respectively. 
The confocal ellipsoidal coordinates are approximately Cartesian near the origin 
and become conical at large radii, i.e., equivalent to spherical coordinates. 
We consider models with a gravitational potential of Stackel form 

a „ „^ _ n>^) EM ^(^) 

»/slA,/.,^J (A-M)(A-z.) (m-^^)(m-A) {u-\){v-tJi)' 

where F{t) is an arbitrary smooth function. This potential is the most general 
form for which the Hamilton-.Jacobi equation separates ([15], [18]) All orbits 
have three exact isolating integrals of motion, which are quadratic in the veloc- 
ities (e.g. [6]). There are no irregular orbits, so that Jeans' theorem is strictly 
valid ([14]), and the distribution function / is a function of the three integrals. 
Therefore, out of the six symmetric second-order velocity moments, defined as 

{viVj){x) = ^jjjviVjf{x,v)d\, (i,i = 1,2,3), (3) 
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with density g, the three mixed moments vanish, and we are left with (vj), {v'^) 
and {vD, related by three Jeans equations. These were first derived by Lynden- 
Bell ([13]), and can be written in the following form ([23]) 

5i(A,/i,z/), (4a) 

g2{X,H,v), (4b) 

gi{\,li,v), (4c) 





Sfifi. 




dX 


2(A-M) 


2{X-iy) 






Sxx 


djj. 




2{^i-x) 




S\x 




du 


2(z/-A) 





where we have defined the diagonal components of the stress tensor 



Srr{X,H,iy) = ^J{X-n){X-v){^l-v)Q{vl) , T = A, /i, J/, (5) 

and the functions gi , g2 and gs depend on the density and potential (2) as 

dVs 



9i{X,n,v) = -^J {X- ij){X-v){ii-v) Q — , (6) 

where 92 and 173 follow from gi by cyclic permutation A — > /x ^ 1/ ^ A. Similarly, 
the three Jeans equations follow from each other by cyclic permutation. The 
stress components have to satisfy the following continuity conditions 

'S'aa(-Q;, -a, v) = Si,^{-a, -a, v) , ^^^(A, -/3) = ^^^(A, -/?) , (7) 

at the focal ellipse (A = /i = —a) and focal hyperbola (p = u = — /?), respectively. 

We prefer the form (5) for the stresses instead of the more common defini- 
tion without the square root, since it results in more convenient and compact 
expressions. In self-consistent models, the density g equals gs, with gs related 
to Vs by Poisson's equation. The Jeans equations, however, do not require self- 
consistency, so that we make no assumptions on the form of g other than that 
it is triaxial, i.e., a function of (A, /i, f), and that it tends to zero at infinity. 



3 The two-dimensional case 

When two or all three of the constants a, /? or 7 in (1) are equal, the triaxial 
Stackel models reduce to limiting cases with more symmetry and thus with fewer 
degrees of freedom. Solving the Jeans equations for oblate, prolate, elliptic disc 
and scale-free models reduces to the same two-dimensional problem ([10], [11], 
[23]), of which the simplest form is the pair of Jeans equations for Stackel discs 



dSxx S^^ 



dX 2{X-n) 
dS^ij, _ Sxx 
djj, 2{ij,-X) 



= .gi(A,/i), (8a) 
= 52(A,m), (8b) 
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with at the foci {X = = —a) the continuity condition 

S\\{-a, -a) = S^i^{-a, -a) . (9) 
In this case the stress components and the functions gi and 52 are 

SrriX.il) ^ y/{X^^ q{vI) (t = A,^), 51 (A, ^) = - ^(V^ , (10) 

where 52 follows from gi by interchanging A <-» /x, and q denotes a surface density. 

The two Jeans equations (8) can be recast into a single second-order partial 
differential equation in cither S\\ or Si^n, which can be solved by employing 
standard techniques like Riemann's method ([5], [23]). However, these standard 
techniques do not carry over to the triaxial case, and we therefore introduce an 
alternative method, based on the superposition of singular solutions. 

We consider a simpler form of (8) by substituting for gi and 52, respectively 
^1=0 and 52 = 5{Xq — X)5{^q—^). We refer to solutions of these simplified 
Jeans equations as singular solutions. Singular solutions can be interpreted as 
contributions to the stresses at a fixed field point (A, /u) due to a source point in 
(Aq, Mo) (Fig- !)■ The full stress at the field point can be obtained by adding all 
source point contributions, each with a weight that depends on the local density 
and potential. Once we know the singular solutions, we can use the superposition 
principle to construct the the solution of the full Jeans equations (8). 

Since the derivative of a step- function TL is equal to a delta-function, it follows 
that the singular solutions must have the form 

5aa = A(A,M)7^(Ao-A)H(A^o-M) , 

(11) 

s^^ = B{x,fM)n{Xo-x)n{fio-fi)-6{Xo-x)n{fio-f^) ■ 

where the functions A and B must solve the homogeneous Jeans equations, i.e., 
(8) with zero right-hand side, and satisfy the following boundary conditions 

-4(Ao,m)= „,/ S(A,Mo) = 0. (12) 
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We solve this two-dimensional homogeneous boundary problem by superposition 

of particular solutions. Wc first derive a particular solution of the homogeneous 
Jeans equations with a free parameter z, which we assume to be complex. We 
then construct a linear combination of these particular solutions by integrating 
over z. We choose the integration contours in the complex z-planc, such that 
the boundary conditions (12) are satisfied simultaneously. The resulting homo- 
geneous solutions are complex contour integrals, which can be evaluated in terms 
of the complete elliptic integral of the second kind, E(w) = \^ y \ — w sin ^, 
and its derivative E'(w), as 



A = 



E(w) 
7r(Ao-/x) ' 



B = - 



2wE'{w) Aq-A QiQ-M 

— , with w = j- — . (13) 

7r(Ao-A) (Ao-Ato)(A-At) 



Wc obtain a similar system of simplified Jeans equations by interchanging the 
expressions for gi and §2 ■ The singular solutions of this simplified system follow 
from (11) by interchanging A /U and Aq Ho at the same time. 

To find the solution to the full Jeans equations (8) at (A,/i), we multiply 
the latter singular solutions and (11) by (7i(Ao,^o) and (72(Ao,/io) respectively, 
and integrate over D = {(Ao,/^o): A < Aq < cxd, ^ < /xq < —a} (Fig. 1). This 
gives the first two integrals of the two equations (14a) and (14b) below. The 
remaining terms are due to the non- vanishing stress at the boundary ^ = —a, 
and are found by multiplying the singular solutions (11), evaluated at ^0 = 
by —S^^{Xo, —a) and integrating over Ag in D. The final result for the solution 
of the Jeans equations (8) for Stackel discs, after using the evaluations (13), is 



oc —a 



Sxx (A, = J j d/^o 



,2wE'{w) , E{w) 

-5i(Ao,Mo)-7 r+52(Ao,Mo) 



A At 



'K{iJ,o-fJ') 7r(Ao-/u) 

E{w) 



- jdXo gi{Xo, m) - ydAo <S',i^(Ao, -a) 



A 

00 —a 



A ^l 



-ff 1 ( Ao , Mo ) — 7T 7 - ff2 ( Ao , Mo ) 



7r(Ao-At) 

2wE'{w) 



, (14a) 

Ho = -oc 



- j dfio g2{X, Mo) + S'^^(A, -a) - JdXo ^'^^^(Ao, -a) 



7r(A-Mo) ' 7r(Ao-A) 

2wE'{w) 
"7r(Ao-A) 



• (14b) 

/io=-a; 



The solution depends on g and Vs through gi and 52- This means that, for given 
Q and Vs, the solution is uniquely determined once we have prescribed S'^^ at 
the boundary ^ = —a. At this boundary, Sw is related to by the first Jeans 
equation (8a), evaluated at m = —C(j up to an integration constant, which is 
fixed by the continuity condition (9). We are thus free to specify either of the 
two stress components at fi = —a. 
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4 The general case 

The singular solution method introduced in the previous section can be extended 
to three dimensions to solve the Jeans equations (4) for triaxial Stackel models. 
Although the calculations are more complex for a triaxial model, the stepwise 
solution method is similar to that in two dimensions. 

We simplify the Jeans equations (4) by setting two of the three functions 
gi, g2 and gs to zero and the remaining equal to 5{\o — X)S{fj,o— n)5{h'o — h'). In 
this way, we obtain three similar simplified systems {i = 1,2,3), each with three 
singular solutions 5'j^'^(A, /x, j^; Aq, ^o, t'o) (t = A,^, t^), that describe the stress 
components at a fixed field point (A, /x, due to a source point in (Aq, ^^o)- 

The singular solutions have a form that is similar to that in the two-dimensio- 
nal case (11). They consist of combinations of step-functions and delta- functions 
multiplied by functions that are the solutions of homogeneous boundary prob- 
lems. The functions that must solve a two-dimensional homogeneous boundary 
problem can be found as in §3, and can be expressed in terms of complete el- 
liptic integrals, cf. (13). The singular solutions in the general case also contain 
three functions A, B and C that must solve the triaxial homogeneous Jeans 
equations, i.e., (4) with zero right-hand side, and satisfy three boundary condi- 
tions. This three-dimensional homogeneous boundary problem can be solved by 
integrating a iwo-parameter particular solution over both its complex param- 
eters, and choosing the combination of contours such that the three boundary 
conditions are satisfied simultaneously. The resulting homogeneous solutions A, 
B and C are products of complex contour integrals, and can be evaluated as 
sums of products of complete /i^/perelliptic integrals. 

To find the solution of the full Jeans equations (4) , we multiply each singular 
solution Sl'^ by ^^(Ao, /io> z^o), so that the contribution from the source point 
naturally depends on the local density and potential. Then, for each coordinate 
T — X,fj,,v, we add the three weighted singular solutions, and integrate over a 
finite volume within the physical region —7 < u < —(3 < /x < — a < A. This 
results in the following general solution of the Jeans equations (4) for triaxial 
Stackel models 




(15) 
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with r = A, jU, v. Whereas the integration hmits A, /u and v are fixed due to the 

position of the field point, the hmits Ag, /ie and v,, arc not, and may be any 
value in the corresponding physical ranges, i.e., Ae G [— a,(X)], € [— /3, — a] 
and Vf. e [—7,-/3], but Ae 7^ —a. The latter choice would lead to solutions 
which generally have the incorrect radial fall-off, and hence arc non-physical, ff 
we choose Ae = 00, there is no contribution from the second line in (15) due to 
vanishing stress at large distance. If we furthermore take = —a and = — 
the integration volume becomes the three-dimensional extension of D (Fig. 1). 

Whereas the volume integral in (15) already solves the inhomogeneous Jeans 
equations (4) , the three area integrals are needed to obtain the correct values at 
the boundary surfaces A = A, , fi = /ie and v = i>e. On each of these siirfaccs the 
three stress components are related by two of the three Jeans equations (4) and 
the continuity conditions (7). Since the (weight) functions gi are known for given 
g and Vs. this means that the solution (15) yields all three stresses everywhere 
in the triaxial model, once one of the stress components is prescribed on the 
three boundary surfaces. If we take Ae = 00 and He = Ve = —P, the contributing 
boundary surfaces reduce to the single {x, z)-plane, containing the long and the 
short axis of the galaxy. This compares well with Schwarzschild ([17]), who used 
the same plane to start his numerically calculated orbits from. 

5 Discussion and conclusions 

Eddington ([9]) showed that the velocity ellipsoid in a triaxial galaxy with a sepa- 
rable potential of Stackel form is everywhere aligned with the confocal ellipsoidal 
coordinate system in which the equations of motion separate. Lynden-Bell ([13]) 
derived the three Jeans equations which relate the three principal stresses to the 
potential and the density. Solutions were found for the various two-dimensional 
limiting cases, but with methods that do not carry over to the general case, which 
remained unsolved. We have presented an alternative solution method, based on 
the superposition of singular solutions (see [23] for details). This approach, un- 
like the standard techniques, can be generalised to solve the three-dimensional 
system. The resulting solutions contain complete (hyper) elliptic integrals, which 
can be evaluated in a straightforward way. 

The general Jeans solution is not unique, but requires specification of prin- 
cipal stresses at certain boundary surfaces, given a separable triaxial potential 
and a triaxial density distribution (not necessarily the one that generates the 
potential). These boundary surfaces can be taken to be the plane containing the 
long and the short axis of the galaxy, and, more specifically, the part that is 
crossed by all three families of tube orbits and the box orbits. 

The set of all Jeans solutions (15) contains all the stresses that are associated 
with the physical distribution functions / > 0, but, as in the case of spherical 
and axisymmetric models, also contains solutions which are unphysical, e.g., 
those associated with distribution functions that are negative in some parts of 
phase space. The many examples of the use of spherical and axisymmetric Jeans 
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models in the literature suggest nevertheless that the Jeans solutions can be of 

significant use. 

While triaxial models with a separable potential do not provide an adequate 
description of the nuclei of galaxies with cusped luminosity profiles and a mas- 
sive central black hole ([7]), they do catch much of the orbital structure at larger 
radii, and in some cases even provide a good approximation of the galaxy po- 
tential. The solutions for the mean streaming motions, i.e., the first velocity 
moments of the distribution function, are helpful in understanding the variety of 
observed velocity fields in giant elliptical galaxies and constraining their intrin- 
sic shapes (e.g. [1], [21], [22]). We expect that the projected velocity dispersion 
fields that can be derived from our Jeans solutions will be similarly useful, and, 
in particular, that they can be used to establish which combinations of viewing 
directions and intrinsic axis ratios are firmly ruled out by the observations. 

It is remarkable that the entire Jeans solution can be written down by means 
of classical methods. This suggests that similar solutions can be found for the 
higher dimensional analogues of (4) , most likely involving hyperelliptic integrals 
of higher order. It is also likely that the higher-order velocity moments for the 
separable triaxial models can be found by similar analytic means, but the effort 
required may become prohibitive. 
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